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ABSTRACT 

Assembly of messenger ribonucleoparticles 
(mRNPs) is a pivotal step in gene expression, but 
only a few molecular mechanisms contributing to 
its regulation have been described. Here, through 
a comprehensive proteomic survey of mRNP 
assembly, we demonstrate that the SUMO pathway 
specifically controls the association of the THO 
complex with mRNPs. We further show that the 
THO complex, a key player in the interplay 
between gene expression, mRNA export and 
genetic stability, is sumoylated on its Hpr1 subunit 
and that this modification regulates its association 
with mRNPs. Altered recruitment of the THO 
complex onto mRNPs in sumoylation-defective 
mutants does not affect bulk mRNA export or 
genetic stability, but impairs the expression of 
acidic stress-induced genes and, consistently, com- 
promises viability in acidic stress conditions. 
Importantly, inactivation of the nuclear exosome 
suppresses the phenotypes of the hpr1 
non-sumoylatable mutant, showing that SUMO-de- 
pendent mRNP assembly is critical to allow a 
specific subset of mRNPs to escape degradation. 
This article thus provides the first example of a 
SUMO-dependent mRNP-assembly event allowing 
a refined tuning of gene expression, in particular 
under specific stress conditions. 



INTRODUCTION 

mRNA packaging into messenger ribonucleoparticles 
(mRNPs) has emerged as a complex and highly regulated 
process which differentially marks processed molecules, 
directs them towards export, and dictates their stability, 
subcellular localization and translation (1). Nuclear 
mRNP assembly initiates with the cotranscriptional 
recruitment of mRNA processing factors and several 
other mRNA-associated proteins, some of which are 
acting as adaptors for the mRNA-export machinery (2). 
In budding yeast, Yral, Nab2 and Npl3 serve as adaptors 
for the Mex67-Mtr2 export receptor, which ultimately 
brings mRNPs to the nuclear pore complex (NPC) 
through interaction with nucleoporins (3). In addition, 
cotranscriptional mRNP assembly has been shown to 
involve THO, which is a conserved tetrameric complex 
comprised of four subunits in yeast — Tho2, Hprl, Mftl 
and Thp2 (4). A pentameric version of THO, including the 
conserved Texl protein, was also recently described (5). 
The THO complex is associated with active genes and 
contributes to the recruitment of several mRNA- 
associated factors, including Sub2, Yral (6,7) and 
SR-like proteins (Gbp2 and Hrbl) (8,9) which together 
with THO form the so-called TREX (TRanscription and 
EXport) complex, and Mex67 (10). Most likely as a con- 
sequence of improper assembly, mRNPs formed in the 
absence of the THO complex are retained in the nucleus 
(7,11,12), and in some cases, degraded by the 
exosome-associated 3—5' exonuclease Rrp6 (13,14). 
THO complex inactivation also alters transcriptional 
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elongation, in particular through G/C-rich, long and/or 
repeat-containing genes (15-17) and provokes genetic 
instability, as revealed by hyper-recombination pheno- 
types (12,18). These phenotypes have been proposed to 
involve R-loops, which are RNA::DNA hybrids that 
accumulate in THO mutants (19,20). 

As highlighted by several reports, changes in mRNP 
composition occur precisely at distinct steps of mRNA 
metabolism (1) whereas improperly assembled mRNPs 
are carefully evicted in the cells (21,22). However, only a 
limited number of spatio-temporal regulations targeting 
mRNP composition have been described. RNA helicase- 
dependent remodeling is best exemplified by the action of 
the DEAD-box ATPase Dbp5, which terminates the 
mRNA export process at the cytoplasmic face of NPCs 
through displacement of mRNA-associated export factors 
such as Mex67 and Nab2 (23). In addition, changes in 
mRNP composition can be driven by mutually exclusive 
interactions involving mRNA-associated proteins and/or 
the mRNA itself: Mex67 recruitment by Yral disrupts the 
Yral-Sub2 interaction in vitro (24), while in mammalian 
cells, mRNA is handed over from adaptors to the export 
receptor upon formation of a ternary complex (25). 
Finally, a number of post-translational modifications tar- 
geting mRNA-associated proteins were also demonstrated 
to modulate mRNP content, impacting on mRNP stabil- 
ity and/or nuclear export (26). Among known examples, 
ubiquitination of Yral disrupts its interaction with 
mRNPs prior to export (27) while Mex67 recruitment 
onto transcribed genes is facilitated by ubiquitination of 
the Hprl subunit of the THO complex (10). 

Protein modification by the small ubiquitin-like 
modifier SUMO, which has emerged as a key regulator 
of multiple cellular processes (28), is also susceptible to 
regulate mRNP assembly. This covalent, reversible modi- 
fication has already been shown to finely tune several 
aspects of mRNA metabolism, including transcription 
and pre-mRNA processing (29,30). In addition, SUMO 
conjugation often leads to conformational changes or 
macromolecular complexes rearrangements (31), and 
could therefore contribute to mRNP remodeling. 
Finally, a number of enzymes of the SUMO pathway 
are associated with nuclear pores, which represent the 
ultimate route for exported mRNPs (30). In particular, 
mRNPs committed for export are docked at the nuclear 
basket of NPCs by virtue of their interaction with myosin- 
like proteins (Mlp) 1 and 2 (32,33), which form a protein 
platform also contributing to the anchoring of the 
conserved SUMO-protease Ulpl (34). Ulpl (SENP2 in 
mammals) is an essential SUMO-isopeptidase, which cata- 
lyzes both the processing of neosynthesized SUMO prior 
to conjugation, and the deconjugation of SUMO from its 
targets (35). A few Ulpl targets have been identified, but 
to date, none of them has been associated with mRNA 
metabolism (35). 

Here, we have investigated the involvement of the 
SUMO pathway, and in particular of the 
SUMO-protease Ulpl, in mRNP assembly. Through a 
comprehensive proteomic analysis of mRNP composition, 
we demonstrate that Ulpl specifically controls the associ- 
ation of the THO complex with mRNPs. We further show 



that the THO complex is sumoylated on its Hprl subunit 
in a Ulpl -dependent manner, and that this modification 
regulates its association with mRNPs. Cells impaired in 
SUMO-dependent THO complex recruitment exhibit 
defective expression of acidic stress-inducible genes and, 
consequently, fail to grow in acidic stress conditions. 
Importantly, exosome inactivation restores both stress- 
induced mRNA levels and viability on stress conditions, 
indicating that defective expression of these mRNAs is 
caused by mRNA degradation triggered by improper 
mRNP assembly. This article thus shows that sumoylation 
of the THO complex regulates mRNP biogenesis, thereby 
contributing to optimal mRNA expression, as revealed 
upon specific stress conditions. 

MATERIALS AND METHODS 

Yeast strains and plasmids 

Unless otherwise indicated, all the strains used in this 
study (listed in Supplementary Table S3) are isogenic to 
S288C and were grown in liquid culture at 25°C and 
shifted for 2h at 30° C. For GAL gene induction, 2% gal- 
actose was added for 2 h to cells grown in glycerol-lactate 
(0.17% YNB, 0.5% ammonium sulfate, 0.05% glucose, 
2% lactate and 2% glycerol) supplemented with the 
required nutrients. For weak acid induction, cells were 
grown on MM4 (minimal medium, pH4) supplemented 
with 60 mM acetic acid as described (36). When indicated, 
cycloheximide (lOOug/ml, Sigma) or methylmethane sul- 
fonate (MMS, 0.02%, Sigma) was added to the medium. 
Construction of plasmids (listed in Supplementary 
Table S4) was performed using standard molecular 
cloning techniques. Sequences of oligonucleotides used 
in this study are available upon request. 

mRNP complexes purification 

mRNP purification from cells expressing proteinA-tagged 
baits was performed essentially as described (37). Briefly, 
frozen grindates were homogenized in nine volumes of 
extraction buffer (20 mM Hepes pH 7.5, HOmM KOAc, 
2mM MgCl 2 , 0.1% Tween-20, 0,5% Triton X-100, 1 mM 
DTT, IX protease inhibitors cocktail, complete 
EDTA-free, Roche and antifoam B, Sigma, 1:5000). 
When indicated, a 10-min treatment with RNAse 
A (Sigma, 150ug/ml) was included. The resulting extract 
was clarified by filtration through 1.6 urn GD/X Glass 
Microfiber syringe filters (25mm, Whatman) and further 
incubated for 30min at 4°C with IgG-conjugated 
magnetic beads. Beads were then washed three times 
with extraction buffer and once with 0.1 M NH 4 OAc, 
0.1 mM MgCl 2 , 0.02% Tween-20. Bound-complexes 
were eluted with 0.5M NH 4 OH, 0.5 mM EDTA, 
lyophilized and resuspended either in SDS-sample buffer 
for SDS-PAGE or in 25 mM ammonium carbonate for 
mass spectrometry analysis. 

Small-scale immunoprecipitations were performed 
using IgG-conjugated magnetic beads or anti-HA 
agarose beads (Pierce). Cells were lysed in the same 
buffer as above by bead beating using a Fastprep 
(Qbiogene). After a 10000-g centrifugation at 4°C for 
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5 min, the soluble extract was incubated with the beads for 
30min (IgG immunoprecipitation) or 90 min (anti-HA 
immunoprecipitation) at 4°C. Beads were washed three 
times with extraction buffer and eluted with SDS sample 
buffer. 

Mass spectrometry analysis 

Digestion was performed overnight at 37°C in the 
presence of 12.5 ug/ml of sequencing grade trypsin 
(Promega, Madison, Wi, USA). Digests were analyzed 
by a LTQ Velos Orbitrap (Thermo Fisher Scientific, San 
Jose, CA) coupled to an Easy nano-LC Proxeon system 
(Thermo Fisher Scientific, San Jose, CA). 
Chromatographic separation of peptides was performed 
with the following parameters: column Easy Column 
Proxeon C18 (10cm, 75 um i.d., 120A), 300nl/min flow, 
gradient rising from 95% solvent A (water — 0.1% formic 
acid) to 25% B (100% acetonitrile, 0.1% formic acid) in 
20 min, then to 45% B in 40 min and finally to 80% B in 
10 min. Peptides were analyzed in the Orbitrap in full ion 
scan mode at a resolution of 30 000 and a mass range of 
400-1800 m/z. Fragments were obtained with a collision- 
induced dissociation (CID) activation with a collisional 
energy of 40%, an activation Q of 0.25 for 10ms and 
analyzed in the LTQ. MS/MS data were acquired in a 
data-dependent mode in which 20 most intense precursor 
ions were isolated, with a dynamic exclusion of 20 s and an 
exclusion mass width of 10 ppm. Data were processed with 
Proteome Discoverer 1.4 software (Thermo Fisher scien- 
tific, San Jose, CA) coupled to an in-house Mascot search 
server (Matrix Science, Boston, MA; version 2.4.1). The 
mass tolerance of fragment ions was set to 7 ppm for pre- 
cursor ions and 0.5 Da for fragments. The following modi- 
fications were considered in mass calculation: oxidation 
(M), phosphorylations (STY). The maximum number of 
missed cleavages was limited to two for trypsin digestion. 
MS/MS data were searched against the SwissProt 
database (11/12/13 release) with the Saccharomyces 
cerevisiae taxonomy. Percolator node was used for 
peptide false discovery rate determination. Peptides were 
validated when their g-value, which represents the 
minimal false discovery rate at which the identification is 
considered correct, was below 1%. 

Chromatin and RNA immunoprecipitation 

Chromatin and RNA immunoprecipitations (ChIP and 
RIP) were performed using a modified procedure (38). 
Briefly, cells were cross-linked in the presence of 1 % for- 
maldehyde for 10 min at 25°C. After quenching with 
100 mM glycine and washing with Tris-buffered saline, 
the cell pellets were resuspended in 1 ml lysis buffer 
(50 mM Hepes pH 7.5, 140 mM NaCl, 1 mM EDTA, 1% 
triton X-100, 0.1% deoxycholate, IX protease inhibitors 
cocktail, complete EDTA-free, Roche) and lysis was 
achieved by bead beating. Cell lysates were collected and 
split in two parts, respectively, used for ChIP or RIP. 
When used for ChIP, the cell lysate was further sonicated 
to shear DNA using a Bioruptor (Diagenode) for two runs 
of 10 min each and solubilized chromatin was retrieved 
through a lOmin-centrifugation at 2500 g to prevent 



differential chromatin fractionation (DCF) (38). When 
used for RIP, the cell lysate was directly spun at 10 000 g 
for 5 min at 4°C. Specific antibodies - anti-Hprl (39); 
anti-RNAP II largest subunit (8WG16, Covance) — were 
added to the relevant soluble extracts and immunopre- 
cipitation was performed by overnight rotation at 4°C. 
Protein-G sepharose beads (GE Healthcare) were added 
to the samples for two additional hours. Washes were as 
follow: twice with lysis buffer, twice with lysis buffer con- 
taining 360 mM NaCl; twice with 10 mM Tris pFt 8, 
250 mM LiCl, 0.5% Nonidet-P40, 0.5% deoxycholate, 
ImM EDTA and once with lOmM Tris-HCl pH 8, 
1 mM EDTA. 

Elution was achieved through a 20-min incubation at 
65°C in the presence of 50 mM Tris pH 8, 10 mM EDTA, 
1 % SDS. The eluate was deproteinized with proteinase K 
(Sigma, 0.2mg/mL) and uncrosslinked for 30 min at 65°C. 
Immunoprecipitated nucleic acids (DNA or RNA) were 
purified with the Qiaquick PCR purification kit (Qiagen) 
and the Nucleospin RNAII kit (Macherey Nagel), respect- 
ively. RNAs were reverse-transcribed with AMV reverse 
transcriptase (Finnzymes). DNA and cDNA were further 
quantified by real-time PCR with a LightCycler 480 
system (Roche) according to the manufacturer's instruc- 
tions. Controls without reverse transcriptase allowed 
estimating the lack of contaminating DNA. RIP (or 
ChIP) values are the ratios between immunoprecipitated 
RNA (or DNA) and input RNA (or DNA) amounts. For 
Hprl ChIP, values were further normalized to the ChIP 
value of RNAP II as described (40). For Hprl RIP and 
ChIP, values were set to 1 for hprl A cells (negative control 
for anti-Hprl immunoprecipitation). 

DCF analyses were performed as described (41). 

SUMO-conjugates isolation 

SUMO-conjugates were isolated from yeast cells express- 
ing a His-tagged version of SUMO under the control of 
the copper-inducible CUP1 promoter. Briefly, 100 OD600 
of cells induced during 4h in the presence of 0.1 mM 
CuS04 were harvested and further processed for nickel 
agarose denaturing chromatography using reported 
methods (42) with the following changes. Cells were 
directly lysed by bead beating in 6M guanidine HC1, 
lOOmM sodium phosphate pH8, lOmM Tris HC1, 0.1% 
TritonX-100, lOmM beta-mercaptoethanol, 50 mM N- 
ethylmaleimide (Sigma). Clarified lysates were incubated 
with Ni-NTA agarose beads (Qiagen) for 2 h at room tem- 
perature. Beads were washed twice with lysis buffer and 
three times with 8M urea, 100 mM sodium phosphate, 
lOmM Tris HC1 pH 6.3 before proceeding to elution. 

Western-blot analysis 

Total protein extraction from yeast cells was performed by 
the NaOH-TCA lysis method (42). Samples were 
separated on 4-12% SDS-PAGE gel and transferred to 
PVDF membranes. Western-blot analysis was performed 
using the following antibodies: rabbit IgG-HRP poly- 
clonal antibody (Dakocytomation), 1:5000; monoclonal 
anti-GFP (Roche Diagnostics), 1:500; monoclonal 
anti-HA (Covance), 1:1000; monoclonal anti-Pabl 
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(Santa-Cruz), 1:1000; polyclonal anti-Yral (6), 1:1000; 
polyclonal anti-Nab2 (32), 1:100000; polyclonal anti- 
Hprl (39), 1:1000; and polyclonal anti-Mex67 (39), 
1:50000. Quantification of signals was performed based 
on serial dilutions of reference samples using the ImageJ 
software. 

Gene expression analysis 

Total RNAs were extracted from yeast cultures using 
Nucleospin RNA II (Macherey Nagel). For RT-qPCR, 
RNAs were reverse-transcribed with AMV reverse tran- 
scriptase (Finnzymes) and cDNA quantification was 
achieved through quantitative real-time PCR as above. 

For microarray analysis, 1 \ig of total RNA was reverse 
transcribed and labeled with Cy5 or Cy3 using the 
amino-allyl protocol described at transcriptome.ens.fr/ 
sgdb/protocols/. A set of 52 519 probes matching the 
ORFs of S. cerevisiae, with 1-15 specific probes per 
gene, was designed using the Teolen software, with 
default parameters (43). This probe design was deposited 
on the earray database (AMADID: 027945) and was used 
to produce Agilent arrays in a 8 x 60 k format. These 
Agilent arrays were hybridized with the labeled cDNA 
using the Agilent hybridization protocol. After washing 
following the supplier's recommendations, the slides 
were scanned using a 2-micron Agilent microarray 
scanner. The images were analyzed using the feature ex- 
traction software and normalized by global lowess using 
Goulphar (44). For each gene, the Cy5/Cy3 ratios corres- 
ponding to the different probes were averaged. The hprl- 
KR versus wt comparison was performed twice using in- 
dependent samples and dye swap. The averaged log 2 of the 
mutant/wild-type ratios and the standard deviation 
between the two replicates were calculated for each gene. 
The genes showing a standard deviation >0.5 were 
removed from the data. This filter eliminated ~ 1 50 
genes which were highly variable between replicates. The 
complete microarray data are available in Supplementary 
Table S2 and in the ArrayExpress database (www.ebi.ac. 
uk/arrayexpress) under accession number E-MTAB-2138. 

Cells carrying LacZ reporters were lysed and assayed 
for f3-galactosidase activity using the NovaBright™ 
(3-galactosidase Chemiluminescent Detection System 
(Invitrogen) following the manufacturer's instructions. 
Relative light units were normalized to the amount of 
cells as estimated by measurement of the OD at 600 nm. 

Miscellaneous 

Fluorescence in situ hybridization (FISH) was carried out 
on fixed cells using a Cy3-(dT)50 probe (45). Wide-field 
fluorescence images were acquired using a DM6000B 
Leica microscope with a lOOx, NA 1.4 (HCX Plan-Apo) 
oil immersion objective and a CCD camera (CoolSNAP 
HQ; Photometries). For Rad52 foci analysis, z-stacks 
sections of 0.2 um were acquired using a piezo-electric 
motor (LVDT; Physik Instrument) mounted underneath 
the objective lens. Images were scaled equivalently and 
3D-projected using MetaMorph 5 (Universal Imaging), 
and further processed with Photoshop CS2 9.0 software 
(Adobe). 



Hyper-recombination rates were determined as the 
frequency of deletions of the chromosomal leu2-k:: 
ADE2-URA3::leu2-k system based on two 2.16-kb leu2 
repeats as previously described (46). 

RESULTS 

A proteomic survey of SUMO function in mRNP 
assembly 

To understand the contribution of sumoylation to mRNP 
assembly and export, we undertook a proteomic analysis 
of mRNPs isolated from cells carrying a thermosensitive 
allele of ULP1 disturbing SUMO-conjugates metabolism 
(ulpl-333, thereafter referred as ulpl ) (47). The ulpl cells 
exhibit both impaired deconjugation and reduced SUMO 
processing at semi-restrictive temperature, thereby 
allowing to assess the global impact of the SUMO 
pathway on mRNP assembly. Using a reported procedure 
maintaining mRNP integrity (37), mRNPs were purified 
from wt or ulpl cells expressing a tagged version of 
the small subunit of the nuclear cap-binding complex 
(Cbp20/Cbc2-protA) and analyzed by SDS-PAGE and 
silver staining (Figure 1A), or by mass spectrometry 
(Supplementary Table S1A). Consistent with previous 
proteomic analysis of Cbp80-associated mRNPs (37), 
Cbc2-associated proteins purified from wt cells included 
a wide range of mRNAs-associated proteins such as 
splicing, 3'-end processing and mRNA export factors 
(« = 85, Supplementary Figure S1A and Table S1A). 
The sensitivity of this proteomic analysis and its selectivity 
for mRNP components were highlighted by the large 
number of recovered interactants (Supplementary 
Table SI A) and the RNAse-sensitive association of most 
identified proteins (Supplementary Table SIB and 
Figure SIB), respectively. These analyses revealed that 
the overall composition of Cbc2-associated mRNPs was 
not drastically altered upon ULP1 loss-of-function, 
indicating that bulk mRNP assembly does not require 
the SUMO pathway (Figure 1A and B and 
Supplementary Table SIB). Remarkably however, out of 
the Cbc2-associated mRNP components whose mass spec- 
trometry scores were strongly altered in the ulpl mutant, 
four corresponded to subunits of the tetrameric THO 
complex (Tho2, Hprl, Mftl, Thp2; Supplementary 
Figure S1A and Table SIB). 

To extend our proteomic survey of mRNP composition, 
we also examined the protein content of nuclear pore- 
associated mRNPs. Proteomic analysis of Mlp2 interact- 
ants unraveled a comprehensive network of proteins 
including all known nucleoporins, spindle pole body com- 
ponents (48), the 26S proteasome (49) and novel partners 
(Figure ID and Supplementary Table S1A). Among them, 
expanding previous findings (32,33), we identified several 
mRNP components (n = 34), whose RNA-dependent 
association with nuclear pores was confirmed by RNAse 
treatment (Supplementary Table SIB and Figure SIC). 
Importantly, the protein contents of Cbc2- and 
Mlp2-associated mRNPs were overlapping but clearly 
distinct (Supplementary Figure S1A and Table S1A), 
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Figure 1. A proteomic survey of SUMO function in mRNP assembly. (A) Cbc2-protA associated mRNPs ('IP') isolated from wt and ulpl cells were 
analyzed by SDS-PAGE followed by silver staining. Position of the Cbc2-protA-tagged protein is indicated. (B) Soluble extracts (left panel, 'inputs') 
and Cbc2-protA-associated mRNPs (mid panel, 'IP') isolated from wt and ulpl cells were analyzed by immunoblotting using the indicated antibodies. 
Decreasing amounts (ul) of a wt reference sample (right panel) were loaded to allow quantification. (C) The relative amounts of Hprl and Mex67 
associated to Cbc2-bound mRNPs (mean ± SD; n = 2) were quantified from (B). Signals recorded in 'no tag' lanes were subtracted and values were 
normalized to the amounts of immunoprecipitated Cbc2 and set to 1 for wt. (D, E and F) same as (A, B and C), for Mlp2-protA-associated proteins. 



thus validating our approach as appropriate to detect 
changes in mRNP composition. In particular, Cbc2- 
associated mRNPs contained multiple splicing factors, 
whereas, consistent with reported nuclear mRNP remodel- 
ing (37), Mlp2-associated mRNPs specifically encom- 
passed late export factors such as Dbp5 or TREX-2 
(transcription and export 2, also referred to as THSC), a 
conserved multi-protein mRNA export factor located to 
the inner face of NPCs (50). The overall composition of 
Mlp2-associated mRNPs did not exhibit strong variations 
upon ULP1 loss-of-function (Figure ID and 
Supplementary Table SIB), indicating that the SUMO 



pathway is not required for mRNP docking at nuclear 
pores. However, mass spectrometry scores of the 
THO complex subunits recovered in the Mlp2-associated 
mRNPs were also decreased upon ULP1 loss-of-function 
(Supplementary Table SIB). This proteomic survey 
thus identified THO complex subunits as a set of 
mRNP-associated proteins whose recruitment to both 
Cbc2- and Mlp2-associated mRNPs could depend on 
Ulpl activity, and therefore, on SUMO metabolism. 
Although Ulpl may also affect the recruitment of other 
mRNA-associated factors (Supplementary Table SIB), we 
focused our subsequent study on the THO complex. 
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SUMO regulates the association of the THO complex 
with mRNPs 

To validate the predictions from our qualitative proteomic 
analysis, we used the same Cbc2- and Mlp2-tagged strains 
for small-scale immunopurifications of mRNPs followed 
by western-blot analysis with a panel of antibodies for 
representative mRNP components. This approach con- 
firmed that ULP1 inactivation decreases the association 
of the Hprl subunit of the THO complex to Cbc2- 
(Figure IB and C) and Mlp2-associated mRNPs (Figure 
IE and F) while other components (Mex67, Yral, Nab2, 
Pabl) remain unaffected. The SUMO pathway is, there- 
fore, specifically required for proper association of the 
THO complex with mRNPs at various stages of the 
export process. 

In parallel, we examined the recruitment of the THO 
complex onto specific target mRNAs using RNA 



immunoprecipitation (RIP) with a specific anti-Hprl 
antibody. Following immunoprecipitation, target mRNAs 
(PMA1, ADH1, YEF3) selected according to previous 
genome- wide analysis of THO complex recruitment (17) 
were quantified by RT-qPCR (Figure 2A). This approach 
detected a uniform association of the THO complex along 
mRNAs that decreased upon inactivation of ULP1 (Figure 
2B), further strengthening the outcome of our mRNPs 
proteomic analysis. 

We next used chromatin immunoprecipitation (ChIP) 
to determine the consequences of ULP1 inactivation on 
the recruitment of the THO complex onto RNAP II- 
transcribed genes (6,7,17). This analysis revealed that the 
association of the THO complex to its aforementioned 
target genes (PMA1, ADH1, YEF3) was not altered by 
ULP1 loss-of-function (Figure 2C). Our data thus indicate 
that this SUMO-protease mainly regulates the association 
of the THO complex with mRNPs. 
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To determine whether these defects were due to 
impaired THO complex assembly, immunoprecipitation 
experiments were carried out in wt and ulpl cells. This 
analysis revealed equivalent coimmunoprecipitation of 
Hprl with Thp2, Mftl or Tho2 subunits, respectively, 
from wt and ulpl cell extracts (Figure 2D). These results 
indicate that the ulpl mutant is not affected for bulk THO 
complex assembly. Taken together, our results therefore 
establish that the SUMO-protease Ulpl, and consequently 
the SUMO pathway, contributes to the recruitment of the 
whole THO complex onto mRNPs. 

The THO complex is sumoylated on the C-terminus of its 
Hprl subunit 

To determine if the THO complex could be a sumoylated 
target of Ulpl, cellular SUMO-conjugates were purified by 
denaturing Ni 2+ chromatography from strains expressing a 
poly-histidine-tagged version of SUMO and HA-tagged 
versions of each of the THO complex subunits. This 
assay did not detect any sumoylated versions of Tho2, 
Mftl and Thp2 subunits (Supplementary Figure S2A). In 
contrary, small amounts of sumoylated Hprl were detect- 
able in SUMO-conjugates fractions from cells coexpressing 
Hprl-HA and His-SUMO (Figure 3A, lane 4, compare 
with lanes 1 and 3). Importantly, sumoylated Hprl was 
also detectable in cells from an unrelated genetic back- 
ground (W303; Supplementary Figure S2A). These 
modified forms were unlikely to correspond to 
ubiquitinated versions of Hprl, which are solely detectable 
in the context of proteasome inhibition (39). In addition, 
these modified forms of Hprl were lost upon simultaneous 
inactivation of the two major SUMO ligases Sizl and Siz2 
(Figure 3B), confirming their sumoylated nature. 
Importantly, the same approach uncovered an accumula- 
tion of sumoylated Hprl in ulpl mutant cells (Figure 3A, 
compare lanes 6 and 4), demonstrating that the 
SUMO-protease Ulpl targets the THO complex for 
desumoylation on its Hprl subunit. 

To further characterize Hprl sumoylation, we went on 
to determine the lysine residues within Hprl that are 
modified by SUMO. For this purpose, we generated 
several plasmid-based Hprl mutants where multiple 
lysines were mutated to arginines to prevent SUMO con- 
jugation without disturbing the charge of the protein 
(Figure 3C) and expressed them in hprlA cells. 
Simultaneous mutation of the first 31 (K1-31R) or the fol- 
lowing 14 (K32-45R) lysine residues did not affect the 
sumoylation level of Hprl (Figure 3D, compare lanes 4, 
5 and 3). In contrast, mutation of the last 20 lysines within 
Hprl C-terminus (K46-65R) completely abolished its 
sumoylation (Figure 3D, compare lanes 6 and 3). Within 
this C-terminal domain, mutation of lysine 46-59 
(K46-50R or K51-59R) did not prevent Hprl sumoylation, 
while simultaneous mutation of lysines 60-65 (K60-65R) 
fully abrogated this modification (Figure 3D, compare 
lanes 10 and 7). 

Similar results were obtained when the hprl-K60-65R- 
HA allele was integrated at the HPR1 locus (Figure 3E). 
Within this domain, the SUMOplot software 
(abgent.com/sumoplot) identified lysine 60 as the unique 



site with a high probability of modification. However, sub- 
stitution of lysines 60 and 61 did not affect Hprl 
sumoylation in our in vivo assay (H. Bretes et al., unpub- 
lished results). This is in agreement with numerous reports 
where mutation of single lysines does not abrogate 
monosumoylation and results in sumoylation at an 
adjacent site (51-53). Taken together, our results thus 
indicate that the THO complex is sumoylated on the 
C-terminus of Hprl. 

THO complex sumoylation occurs independent of its 
ubiquitination 

Ubiquitination was previously demonstrated to target 
Hprl for proteasomal degradation (39). In conditions 
where sumoylated Hprl readily accumulates (e.g. in the 
ulpl mutant, Figure 3A), the degradation rates of Hprl or 
other THO complex subunits were not altered (Figure 4A 
and B). Similarly, in conditions where Hprl sumoylation 
is prevented (e.g. in the sizl siz2 SUMO ligase mutant, or 
in the hprl-K60-65R mutant, see Figure 3B and E), the 
global turnover of Hprl remained unaffected (Figure 4C 
and D). These results thus demonstrate that sumoylation 
of Hprl neither prevents, nor stimulates its ubiqui tin-de- 
pendent degradation. 

In contrast, simultaneous mutations of lysines 51-59 led 
to a strong accumulation of the Hprl protein (Figure 4E 
and F; see also Figure 3D, compare lane 9 and 7, or 19 
and 17), indicating that these lysines likely contribute to 
Hprl ubiquitin-dependent proteasomal degradation. 
Importantly, the Hprl-K51-59R mutant protein was still 
sumoylated (Figure 3D), showing that Hprl 
ubiquitination does not modulate its sumoylation. Of 
note, Hprl C-terminus (residues 548-752) was identified 
as a partner of the ubiquitin-associated domain of Mex67 
in the two-hybrid assay (10) and Hprl degradation was 
shown to be slowed down upon deletion of its last 88 
aminoacids (39). Together with our results, these data 
strongly suggest that within lysine 51-59, lysines 58 and 
59 are the major residues involved in ubiquitin-dependent 
degradation of Hprl (Figure 4G). 

Interestingly, secondary-structure-prediction programs 
detect with a high probability the presence of an intrinsic- 
ally disordered domain in this region of the Hprl protein 
(residues 667-752; Supplementary Figure S2B). Such 
domains have been demonstrated to be particularly 
prone to modification by ubiquitin (54) and ubiquitin- 
like modifiers (55). Together, our data reveal that 
ubiquitination and sumoylation occur independently on 
distinct residues within a short and likely disordered 
domain at the C-terminus of Hprl. 

THO complex sumoylation regulates its association with 
mRNPs 

Although HPR1 deletion differentially impacted viability 
depending on the temperature and the genetic back- 
ground, these growth defects were fully rescued by expres- 
sion of the Hprl-K60-65R protein (Figure 5 A). Since this 
allele does not impair Hprl ubiquitination and comple- 
ments HPR1 deletion, it was further used as a specific 
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Figure 5. THO complex sumoylation regulates its association with mRNPs. (A) Serial dilutions of the indicated strains were grown on YPD plates at 
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indicated strains by RIP using an anti-Hprl antibody followed by RT-qPCR using specific primer pairs. The positions of the primers along the 
PMA1 gene are indicated. RIP values (mean ± SD; = 3) are the ratios between immunoprecipitated and input RNA amounts. Values were set to 1 
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non-sumoylatable hprl mutant in our next studies (subse- 
quently referred to as hprl-KR). 

To evaluate the impact of Hprl sumoylation on THO 
complex association with mRNAs, we analyzed the 
recruitment of the hprl-KR mutant protein into 
Cbc2-associated mRNPs. This revealed a reduced associ- 
ation of the non-sumoylatable hprl-KR protein with 
mRNPs as compared to wt Hprl (Figure 5B, compare 
lanes 5 and 6; Figure 5C). To better quantify this 
decreased recruitment, we performed RIP in the same 
mutants (Figure 5D). This approach confirmed the 
reduced association of the non-sumoylatable hprl-KR 
mutant with its target mRNAs (Figure 5D). Both 
approaches further revealed that unlike in HPR1 wt 
cells, the decreased recruitment observed in hprl-KR 
cells was not further enhanced by ULP1 inactivation 
(Figure 5B, compare lanes 6 and 8; Figure 5C and D). 
This strongly supports that a sumoylation/desumoylation 
cycle occurring on the C-terminal domain of Hprl 
contributes to its recruitment onto mRNPs. 

Importantly, sumoylation did not impair the bulk 
association of Hprl with the THO complex, as revealed 
by the similar coimmunoprecipitation of both Hprl and 
hprl-KR with HA-tagged Thp2, Mftl or Tho2 
(Figure 5E). Although we cannot formally exclude that 
the assembly of the specific mRNP-associated fraction of 
the THO complex could be altered in hprl-KR and ulpl 
mutants, our data indicate that Hprl sumoylation likely 
regulates the recruitment of the entire THO complex onto 
mRNPs. 

THO complex sumoylation does not impact on genetic 
stability and bulk mRNA export 

In view of the described functions of the THO complex, 
we next analyzed the consequences of impaired Hprl 
sumoylation on genetic stability, mRNA export and 
mRNA expression. Unlike THO-inactivation mutants 
(hprl A) (18), mutants affecting THO complex 
sumoylation (hprl-KR, ulpl) do not display increased 
recombination between direct chromosomal repeats 
(Supplementary Figure S3 A). In addition, impairing 
sumoylation of Hprl does not lead to synthetic lethality 
when combined with DNA damage checkpoint mutants 
(Supplementary Figure S3B and H. Bretes et al., unpub- 
lished results), a phenotype previously reported for THO 
null mutants (56). Finally, unlike hprlA cells (56), 
hprl-KR cells do not exhibit an increased occurrence of 
DNA double-strand break-repair foci labeled by Rad52 
(Supplementary Figure S3C). Taken together, these data 
show that THO complex sumoylation does not impact on 
genetic stability. 

Recruitment of known THO complex-dependent 
mRNP components such as Mex67 was not detectably 



affected by changes in Hprl sumoylation in the ulpl or 
hprl-KR mutants (see Figures IB and E and 5B). In agree- 
ment with these biochemical data, the ulpl and hprl-KR 
mutants did not exhibit bulk mRNA export defects, as 
probed by FISH analysis of poly-A+ RNAs 
(Supplementary Figure S4A). These results indicate that 
while THO complex sumoylation regulates its own asso- 
ciation with mRNPs, it does not further affect mRNP 
export. 

THO complex sumoylation is required for the expression 
of stress-inducible genes and for viability upon stress 
conditions 

To gain insight into the impact of THO complex 
sumoylation on gene expression, we first analyzed the 
expression of the bacterial LacZ gene, a reporter typically 
affected in THO mutants. Beta-galactosidase assays 
revealed a mild yet significant reduction of LacZ expres- 
sion in hprl-KR cells as compared to wt (Supplementary 
Figure S4B). 

This result prompted us to characterize the impact of 
THO complex sumoylation on endogenous gene expres- 
sion. Total RNAs were extracted from wt and hprl-KR 
cells and transcriptome profiles were analyzed by micro- 
array. Preventing Hprl sumoylation (hprl-KR) did not 
cause extensive changes in gene expression at the 
genomic level (Figure 6A). Indeed, only 41 genes exhibited 
mRNA levels at least 1.5-fold above or below the wt 
values in hprl-KR cells (Supplementary Table S2), thus 
differing from the reported genome-wide impact of THO 
inactivation (811 genes affected in hprl A cells) (17). In 
addition, only a very limited number of genes exhibited 
consistent changes between two replicate experiments. 
Among them, the only clear functional enrichment was 
observed for genes involved in the adaptive response to 
chemical stresses caused by lowered pH (Haal regulon). 
These genes (YR02, TP02, TDA6, HSP30, YLR297w, 
YDR222w, TP03, indicated on Figure 6A) were therefore 
selected for further studies. 

To confirm that Hprl sumoylation is required for the 
expression of this subset of genes, wt and hprl-KR cells 
were grown in conditions of induction of the Haal 
regulon, e.g. minimal media (MM4) supplemented with 
acetic acid (36). In agreement with previous reports (57), 
acid treatment lead to an increased expression of the 
considered genes in wt cells, as revealed by RT-qPCR 
(YR02, TDA6, TP02, TP03 panels; Figure 6B). 
Importantly, hprl-KR cells specifically exhibited reduced 
mRNA levels for these genes (Figure 6B, compare with the 
ACT1 panel). 

Expression of genes of the Haal regulon is essential for 
cell survival upon acidic stress (36). Strikingly, impairing 
Hprl sumoylation (hprl-KR) similarly affected viability 



Figure 5. Continued 

for hprl A cells (negative control for anti-Hprl immunoprecipitation). n.d., not detectable. */><0.05 (Student's /-test, compared to wt). (E) Soluble 
extracts from wt, THP2-HA, THP2-HA hprl-KR, MFT1-HA, MFT1-HA hprl-KR, TH02-HA or TH02-HA hprl-KR cells were used for anti-HA 
immunoprecipitation. Soluble extracts (left panels, 'inputs') and immunoprecipitated material (right panels, 'a-HA IP') were analyzed by western 
blotting using anti-Hprl (top panels) and anti-HA (lower panels) antibodies. Note that Hprl levels are slightly decreased in the presence of HA tags 
on Mftl or Tho2; however, each of these HA-tagged subunits coimmunoprecipitate equivalent amounts of Hprl from wt and hprl-KR extracts. 
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Figure 6. THO complex sumoylation is required for the expression of stress-inducible genes and for viability upon stress conditions. (A) Transcriptome 
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on acetic acid-containing media (Figure 6C). Notably, the 
growth defect of hprl-KR cells recorded on acid-contain- 
ing media was also observed in an unrelated genetic back- 
ground, albeit to a lower extent (W303; Supplementary 
Figure S5A). However, this phenotype did not reflect a 
general sensitivity of the hprl-KR mutant to stress condi- 
tions (Supplementary Figure S5B). In particular, unlike 
THO inactivation mutants (13), hprl-KR cells did not 
display heat sensitivity (Figure 5A) or 3'-end degradation 
of the heat-shock induced HSP104 transcript 
(Supplementary Figure S5C). Taken together, our results 
thus establish that Hprl sumoylation is critical for the 
optimal expression of a subset of stress-inducible genes 
(e.g. the Haal regulon) and subsequently, for viability in 
the corresponding stress conditions. 

Improperly assembled mRNPs are targeted for 
degradation by the exosome in THO sumoylation mutants 

We first confirmed that sumoylation of Hprl also 
regulates THO complex association with acid-induced 
mRNAs. RIP analyses indeed revealed that both 
hprl-KR and ulpl cells exhibit decreased association of 
the THO complex with TP02 and YR02 mRNAs 
(Figure 6D). 

The expression of this subset of genes is notably 
regulated by the Haal transcriptional activator (58). 
However, hprl-KR cells were not affected for the expres- 
sion of the HAA1 mRNA (Figure 6B). In addition, the 
effect of the hprl-KR mutant was observed both in 
uninduced and acid-induced conditions (Figure 6B). 
Furthermore, RNAP II recruitment along acid-induced 
genes was not affected in hprl-KR cells (Supplementary 
Figure S5D). These results therefore suggest that the 
observed effects on the Haal regulon in hprl-KR cells 
are post-transcriptional. 

To assess whether improper mRNP assembly triggers 
increased degradation of mRNAs of the Haal regulon, 
the hprl-KR allele was combined with the inactivation of 
the Rrp6 subunit of the exosome. RRP6 deletion clearly 
suppressed the mRNA expression defect in hprl-KR cells, 
restoring the amounts of YR02, TDA6, TP02 and TP03 
acid-induced mRNAs to levels comparable with wt 
(Figure 6B). Importantly, while RRP6 inactivation 
rescued TP02 and YR02 mRNAs levels in hprl-KR 
cells, it did not restore Hprl association to these 
mRNAs (Figure 6D), further strengthening the fact that 
exosomal degradation occurs as a consequence of 
decreased Hprl recruitment. Moreover, RRP6 deletion 
suppressed the strong viability defect of hprl-KR cells on 
acid-containing media (Figure 6C). Taken together, our 
data thus establish that in THO sumoylation mutants, 
improperly assembled mRNPs of the Haal regulon are 
targeted for degradation by the exosome. 

In THO-inactivation mutants, exosome-mediated 
degradation of HSP transcripts was previously reported 
to be associated with defective mRNP release from the 
transcription site. This leads to the persistence of the cor- 
responding genomic loci at the vicinity of the NPC, giving 
rise to a detectable chromatin heterogeneity referred to as 
DCF (38). While we detected DCF formation both at 



HSP104 and acid-induced genes loci upon THO inactiva- 
tion, this phenotype was not recorded in hprl-KR cells 
(Supplementary Figure S5E). Loss of THO complex 
sumoylation therefore triggers defects in mRNP 
assembly and stability for acid-induced mRNAs without 
detectable mRNP retention (as probed by DCF). 



DISCUSSION 

While an increasing number of post-translational modifi- 
cations have been reported to influence mRNP assembly 
(26), we describe for the first time the involvement of 
SUMO in this process. Indeed, we show that the Hprl 
subunit of the THO complex is sumoylated on its 
C-terminus, and that its sumoylation/desumoylation regu- 
lates the association of the THO complex with mRNPs. 
Distinct mechanisms may account for the intriguing 
finding that both the non-sumoylatable (hprl-KR) and 
the hyper-sumoylated (ulpl) mutants display decreased 
THO recruitment on mRNPs (Figure 5). 

On one hand, SUMO was already reported to modify 
the interaction properties of its target proteins with nucleic 
acids, either DNA (59) or RNA (52,60). Hprl 
sumoylation could directly interfere with the nucleic-acid 
binding properties of the THO complex, for example by 
preventing the recently demonstrated interaction between 
Tho2 C-terminus and DNA/RNA (5). Interestingly, Sizl 
and Siz2, the SUMO ligases for Hprl (Figure 3B), have 
been reported to bind DNA through their SAP domain 
(61,62). For genes at the nuclear periphery, an Hprl 
sumoylation-desumoylation cycle driven by chromatin- 
bound SUMO ligases and NPC-bound Ulpl could be 
necessary for mRNA recognition by the THO complex. 

On the other hand, sumoylation-desumoylation could 
be part of complex formation, as suggested by the fact that 
several transcription factors are recruited into co-repres- 
sor complexes in a SUMO-dependent manner, but do not 
appear to be sumoylated in the final assembly (63). 
Likewise, Hprl sumoylation could modulate the inter- 
action of the THO complex with an unknown protein 
partner involved in its recruitment onto mRNPs. 
Interestingly, sumoylation of intrinsically disordered 
domains such as the one encompassing Hprl SUMO 
acceptor sites has been proposed to create novel inter- 
action surfaces (51) which, in the case of the THO 
complex, would allow mRNP binding. Following this 
initial recruitment, desumoylation of Hprl could 
remodel the mRNP, resulting in stable association of the 
THO complex. Of note, these different scenarios would 
involve a transient sumoylation of the THO complex, 
which is consistent with the detection of low levels of 
sumoylated Hprl in wt cells (Figure 3A). 

Our data further revealed that sumoylation of Hprl is 
largely independent from its ubiquitination-dependent 
turnover (Figure 4). However, other post-translational 
modifications could also control THO complex activity 
as reported for other mRNA-associated factors (26). For 
instance, DNA damage-induced phosphorylation of the 
THO complex was reported to trigger its dissociation 
from RNA in mammalian cells (64). Spatio-temporal 
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coordination of these post-translational modifications is 
also expected to be critical for mRNP fate. Strikingly, 
the only conserved enzymes reported to associate with 
nuclear pores are the SUMO-protease Ulpl and 
the RNA helicase Dbp5, which both target 
mRNPs-associated factors. Mechanisms ensuring a local 
activation of Dbp5 at the cytoplasmic face of nuclear 
pores have been extensively documented (23). While it is 
tempting to speculate that Ulpl -dependent remodeling 
could occur at the nuclear basket, a free fraction of 
Ulpl may also perform this task in the nucleoplasm. In 
both hypotheses, physiological conditions of disturbed 
Ulpl localization (65,66) could influence mRNP 
composition. 

In agreement with an overall normal mRNP compos- 
ition, ulpl and hprl-KR mutants do not display any 
mRNA export defects (Supplementary Figure S4), sug- 
gesting that the loading of mRNA export factors onto 
mRNPs is fully supported by the chromatin-associated 
THO complex fraction, or compensated by redundant 
pathways (67,68). In contrast, we have demonstrated 
that SUMO-dependent control of THO complex recruit- 
ment onto mRNAs is critical for acidic stress-induced 
gene expression. Indeed, hprl-KR cells exhibit strongly 
reduced levels of mRNAs of the acid-induced, Haal -de- 
pendent regulon and do not survive on acid-containing 
media. Altered expression of this subset of genes thus 
likely accounts for the previously reported identification 
of the hprlA and thp2A THO inactivation mutants in 
several genome-wide screens for weak acid sensitive 
strains (69-71). The fact that expression of the HAA1 
regulator itself and RNAP II recruitment along 
acid-induced genes were not affected in hprl-KR cells sug- 
gested that the observed effects on the Haal regulon could 
be post-transcriptional. Consistently, inactivation of the 
exosome (rrp6A) was shown to restore acid-induced 
mRNAs levels in the hprl-KR mutant. These results 
indicate that in wt cells, SUMO-dependent recruitment 
of the THO complex allows this specific subset of 
mRNPs to escape exosomal degradation, either by 
directly protecting the mRNA, or by contributing to the 
recruitment of other mRNP-associated factors that would 
not have been uncovered in our proteomic analysis of bulk 
mRNP composition. Interestingly, sub-optimal mRNP 
assembly has already been proposed to trigger 
exosome-mediated degradation of heat shock-induced 
HSP mRNAs in THO deletion mutants (13,72). Of note, 
THO null mutants are extremely sick in combination with 
exosome inactivation (13), suggesting that the failure to 
degrade improperly assembled mRNPs is detrimental. In 
contrast, hprl-KR cells, which have more limited defects in 
mRNP assembly, are not sensitive to exosome inactivation 
under normal growth conditions. Moreover, RRP6 
deletion rescues the lethality of hprl-KR cells on acid-con- 
taining media, likely by restoring normal levels of mRNAs 
of the Haal regulon. 

mRNAs which are affected upon impairment of THO 
function share the common feature of being strongly ex- 
pressed, in particular in response to stress (13,17,38,41). 
Indeed, strong bursts of mRNP production might be 
highly sensitive to any alterations in mRNP packaging 



arising upon THO dysfunction. However, loss of Hprl 
sumoylation does not affect the levels of any abundant 
or stress-induced mRNAs. Additional features of acid- 
induced genes of the Haal regulon may render their 
mRNAs more sensitive to degradation in hprl-KR cells. 
This might include particular chromatin environment, 
nuclear positioning or secondary RNA structures that 
could also favor exosomal degradation upon improper 
mRNP packaging. While the initial identification of this 
specific subset of genes in our transcriptome analysis is 
due to the fact that they are already down-regulated in 
uninduced hprl-KR cells, other mRNAs, not expressed 
in our conditions, may also require Hprl sumoylation 
for proper expression. In the future, detailed analysis of 
mRNP composition for specific mRNAs will be required 
to further refine the role of THO complex sumoylation in 
mRNP assembly, both in normal and environmentally 
challenged cells. 
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